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I.  INTRODUCTION 


The  study  of  the  active  solar  corona  requires  information  about 

such  thermodynamic  variables  as  the  plasma  temperature,  density,  and 

pressure  as  a  function  of  position.  X-ray  spectra  allow  the  dlstrlbu- 

tion  of  emission  measure  /  n  dV  as  a  function  of  temperature  to  be 

© 

determined  in  flaring  (McKenzie  and  Landecker  1981)  and  nonflaring 
active  regions.  V  is  the  volume  and  ne  the  electron  density.  The 
usefulness  of  the  emission  measure  is  limited  because  ne  appears  in 
the  second  power.  In  general,  such  quantities  as  the  electron  densl- 
ty,  ne,  or  the  pressure,  nekT,  cannot  be  derived  from  emission  measure 
alone.  A  direct  measurement  of  tx^  is  required. 


A  small  number  of  density-diagnostic  X-ray  line  ratios  are  avail¬ 
able.  Gabriel  and  Jordan  (1969)  showed  that  the  flux  ratio,  R,  of  the 
ls^  -  ls2s  forbidden  line  to  the  ls^  *Sq  -  ls2p  inter- 
combination  line  in  the  helium-like  ions  would  be  density-sensitive 
for  densities  exceeding  a  species-dependent  lower  limit  n£*  (in  this 
paper  we  consider  flux  ratios  where  flux  density  has  units  of  photons 
cm~^  s'* ).  Hie  ratio  in  0  VII  (McKenzie  et  al.  1980a;  Doschek  et  al. 
1981)  has  been  used  to  derive  coronal  densities  in  solar  flares.  R 
for  Ne  IX,  Mg  XI,  and  Si  XIII  may  also  prove  to  be  useful,  but  for 
ions  heavier  than  these  ne*  is  almost  certainly  larger  than  ne  ever 
gets  in  the  corona.  In  addition,  a  line  ratio  in  Fe  XXI  (Mason  et  al. 
1979)  provides  a  density  diagnostic  for  ne  >  10  cm  .  McKenzie  and 


Landecker  (1982)  have  shown  that  the  same  ratio  in  Ca  XV  is  useful  for 
measuring  densities  that  are  commonly  present  in  flares  at  tempera- 


tures  of  ~  4  x  10^  K.  The  Fe  SCI  diagnostic  has  not  yet  been  used 
because  at  ~  10  ^  K,  the  temperature  at  which  Fe  XXI  is  aost  abundant, 
densities  are  apparently  rarely  high  enough  to  yield  a  measurable 
effect  (McKenzie  et  al.  1980b). 

The  0  VII  density  diagnostic  is  especially  Interesting  because  it 
offers  the  capability  of  measuring  ne  in  at  least  some  nonflaring 
active  regions.  The  0  VII  lines  are  produced  at  ~  2  x  106  K,  in  the 
range  of  temperatures  commonly  observed  in  coronal  active  regions. 
Furthermore  ne*(0  VII)  »  3  x  10®  cm~^,  a  density  that  is  sometimes 
exceeded  in  the  nonflaring  corona.  A  major  aim  of  the  present  paper 
is  to  examine  the  usefulness  of  the  0  VII  ratio  in  the  nonflaring 
corona.  To  do  this  we  analyze  data  from  the  SOLEX  B  spectrometer 
aboard  the  U.  S.  Air  Force  P78-1  satellite  (Landecker,  McKenzie,  and 
Rugge  1979;  McKenzie  et  al.  1980b).  The  resolving  power  of  the  spec¬ 
trometer  is  not  high  enough  to  permit  a  similar  study  of  R  for  Ne  IX. 

The  helium-like  line  ratio, 

F(ls2  XS0  -  ls2s  3S1)  +  F(ls2  -  ls2p  3Pj) 

G ■  21  1  * 

FUs*  AS0  -  ls2p  ‘p^ 

where  F  is  flux,  is  useful  as  an  indicator  of  the  state  of  ionization 
equilibrium  during  flares  (Acton  and  Brown  1978;  see  also  Gabriel  and 
Jordan  1969,  note  added  in  proof).  In  addition  the  resonance  line 
(denominator)  is  depleted  by  resonance  scattering  in  the  corona  (Acton 
1978),  so  G  can  be  used  as  a  measure  of  the  coronal  column  density  of 
0  VII.  For  0  VII  the  resonance  line  is  most  strongly  emitted  at  temp- 


erature  T  ~  2  x  106  K,  but  the  population  of  the  absorbing  0  VII  Ions 
Is  highest  at  ~  I  x  10^  K  (Jacobs  et  al.  1978).  Thus  G  can  be  used  as 
a  probe  of  the  quiet  corona  between  the  observer  and  the  emitting 
active  region.  The  spectra  of  both  0  vn  and  Ne  IX  show  resonance 
scattering  effects.  We  examine  G  as  a  plasma  diagnostic  for  both 
species  In  §  III. 


II.  THE  0  VII  DENSITY  DIAGNOSTIC 

The  electron  density  In  a  2  x  10^  K  plasma  can  be  derived  from 
the  ratio  R  for  0  Vll  as  follows  (Pradhan  and  Shull  1981): 


ne  -  3.1  x  1010  (^2.  -  i)  a,,"3 


(2) 


R0  Is  the  value  taken  by  R  In  the  limit  of  low  densities.  It  Is  cus¬ 
tomary  to  assume  that  densities  for  which  RQ/R  -  1  >  0.1  are  measur¬ 
able  (Gabriel  and  Jordan  1969),  but  this  may  not  always  be  true  since 
statistical  errors  In  R  and  uncertainties  in  RQ  both  contribute  to  the 
error  In  nfi.  Differentiating  equation  (2),  we  find 
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3.1  x  10 
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cm 


-3 


(3) 


Even  if  statistical  uncertainties  In  R  are  negligible,  for  active 
region  densities  of  a  few  times  107  cm  ,  a  10%  uncertainty  in  RQ  can 
result  In  very  large  errors  in  the  derived  ne.  Obtaining  an  accurate 
value  for  RQ  is  of  primary  importance  in  applying  the  density  diag¬ 
nostic  to  nonflaring  plasmas. 
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Gabriel  and  Jordan  (1972)  calculated,  for  0  VII,  RQ  -  3.6,  and 
this  value  has  largely  been  accepted  ever  since.  More  recently 
Pradhan  and  his  colleagues  have  recalculated  RQ  taking  into  account 
autoionizing  resonances  in  the  collisional  excitation  cross  sections 
(Pradhan,  Norcross,  and  Hummer  1981;  Pradhan  and  Shull  1981).  Pradhan 
and  Shull  find  RQ(0  VII)  *  3.95.  They  estimate  that  the  overall  error 
in  their  line  ratios  (for  a  number  of  helium-like  ions)  should  not 
exceed  25%,  but  recently  Pradhan  (1981)  has  estimated  that  the  errors 
for  ions  other  than  Fe  XXV  and  Ca  XIX  ("particularly  0  VII”)  should  be 
"well  below  10%.” 

Prior  to  the  P78-1  launch  on  1979  February  24,  a  substantial  data 
set  of  R(0  VII)  measurements  existed.  Rugge  and  Walker  (1970,  1971) 
compiled  data  from  the  OV1-10  satellite  for  90  days  in  1966/1967  and 
from  the  0V1-17  satellite  for  a  few  days  in  1969.  They  found  mean 
values,  R,  of  3.3  and  3.2,  respectively,  in  the  two  data  sets.  The 
data  can  be  criticized  in  that  the  spectral  resolution,  particularly 
in  the  0V1-10  data,  does  not  allow  the  resonance  and  intercombination 
lines  to  be  resolved  well,  and  the  background  in  the  relevant  wave¬ 
length  region  is  consequently  difficult  to  estimate.  Substantial 
systematic  errors  may  be  in  effect.  Acton  et  al.  (1972)  and  Acton  and 
Catura  (1976)  report  measurements  of  R(0  VII)  by  rocket-borne  spectro¬ 
meters  with  much  better  spectral  resolution.  In  each  case  the  ratios 
are  distributed  about  a  mean  of  ~  3.6.  The  Acton-Catura  (1976)  data 
include  one  region  in  which  R/R  ■  1.25.  The  authors  Interpret  this  as 
a  positive  measurement  of  the  density.  Furthermore,  they  point  out 
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Chat,  as  a  result  of  the  co-existence  in  the  spectrometer  field-of- 
view  of  high-  and  low-density  regions,  the  R  value  of  3.60  was  below 
the  low-density  limit  value,  RQ.  One  measurement  by  Acton  ^t_  al. 
(1972)  and  one  by  McKenzie  et  al.  (1978)  are  of  particular  interest  in 
that  active  regions  were  excluded  from  the  respective  fields  of 
view.  The  values  of  R  from  these  measurements  should  be  close  to 
Rq.  Acton  et  al.  (1972)  found  R  ■  3.4  ±0.1  and  McKenzie  et  al. 
(1978)  found  R  ■  4.2  ±  0.2.  The  former  value  is  surprisingly  low;  in 
fact  it  is  lower  than  or  equal  to  R  in  all  six  active  regions  or  com¬ 
plexes  of  active  regions  reported  in  the  same  paper.  Perhaps  unde¬ 
tected  dense  material  was  within  the  field  of  view.  The  McKenzie  et 
al.  (1978)  result  is  consistent  with  the  Pradhan  and  Shull  (1981) 
value  of  RQ  but  not  with  that  of  Gabriel  and  Jordan  (1972). 

In  summary,  theoretical  values  of  RQ  range  from  3.6  to  3.93,  the 
latter  value  arising  from  vwsre  recent  and  sophisticated  calcula¬ 
tions*  Observations  of  nonflaring  plasmas  have  yielded  values  of  R 
from  ~  3.2  to  ~  4.2.  The  best  of  the  rocket  observations  appear  to  be 
those  of  Acton  and  Catura  (1976),  which  combine  high  counting  rates 
with  spectral  resolution  adequate  to  resolve  the  three  0  VII  lines. 
The  McKenzie  et  al.  (1978)  result  comes  close  to  being  a  direct  meas¬ 
urement  of  Rq  and  combines  high  counting  rates  with  adequate  spectral 
resolution  (somewhat  worse  than  that  of  Acton  and  Catura  (1976)). 
Both  the  Acton-Catura  and  McKenzie  et  al.  results  support  R0  >  3.60. 

The  SOLEX  B  data  set  to  be  considered  here  consists  of  29  meas¬ 
urements  of  R  (and  G)  made  between  1979  April  5  and  July  22.  Each 
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measurement  uses  data  from  a  sum  of  ~  10  successive  spectral  scans* 
Individual  spectral  scans  take  either  84  s  or  168  s  to  traverse  the 
Bragg  angle  range  of  17.4°  -  61.7°  In  30". 2  steps.  The  spectrometer 
is  collimated  to  a  60"  (FWHM)  square  field  of  view,  and  in  these  meas¬ 
urements  a  rubidium  acid  phthalate  crystal  (RAP;  2d  ■  26.121  A)  was 
used.  With  possible  exceptions  to  be  discussed  below,  the  spectra  are 
of  nonflaring  active  regions,  but  the  collimator  restricts  the  field 
of  view  so  that  the  brightest  or  densest  parts  of  the  regions  are  not 
always  viewed.  The  0  VII  lines  are  barely  detectable  outside  active 
regions.  Figure  1  shows  two  sample  (slow  scan)  0  VII  summed 
spectra.  The  spectra  are  comparable  to  the  Acton  and  Catura  (1976) 
rocket  measurements  in  both  spectral  resolution  and  counting  statis¬ 
tics.  The  forbidden  line  integrated  count  less  background  is  typical¬ 
ly  several  thousand.  The  SOLEX  spectra  have  the  advantage  of  finer 
colllmatlon  than  the  1.3'  one-dimensional  collimator  used  in  the 
Acton-Catura  observations. 

We  first  attempted  to  fit  the  SOLEX  data  in  the  wavelength  range 
shown  in  Figure  1  to  a  model  consisting  of  the  continuum  and  seven 
spectral  lines:  the  three  0  VII  lines,  the  Ca  XVI  2p  ^P^/2  “  3d  ^>3/2 
line  at  21.444  A  (McKenzie  and  Landecker  1982),  and  three  unidentified 
weak  lines  at  22.03  A,  21.94  A,  and  21.97  A.  The  last  two  lines  men¬ 
tioned  were  marginally  detected  in  flare  spectra;  they  had  a  negli¬ 
gible  effect  on  the  fits.  The  background  was  fitted  using  a  two- 
parameter  fit  which  was  linear  in  the  calculated  spectrometer  sensi¬ 
tivity  as  a  function  of  wavelength.  A  Voigt  line  profile  was  used 
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Figure  1.  X-Ray  Spectra  of  Two  Solar  Active  Regions  Showing  Strong 
Lines  of  0  VII.  The  spectra  were  chosen  to  illustrate 
the  variability  in  the  data.  The  upper  spectrum  has 
R  -  3.4  and  G  -  1.8,  while  the  lower  has  R  •  3.9  and 
G  -  1.2. 


5 


(Finn  and  Muggiest  one  1965)  with  the  line  positions  and  widths  as  free 

2 

parameters.  In  only  ten  of  the  27  spectra  fitted  using  a  minimum— x 

2  2 

criterion  was  the  x  acceptable  at  the  95%  x  probability  level  (with 

2 

442  degrees  of  freedom,  x  <  491  was  acceptable).  Ue  attribute  this 
failure  to  difficulty  in  modeling  the  line  profiles  of  the  strong  0 
VII  lines.  Resonance  scattering  effects,  to  be  discussed  in  §  III, 
can  significantly  distort  the  resonance  line  profile,  especially  near 
line  center. 

Since  the  results  from  the  spectral  fits  were  unsatisfactory,  we 
used  another  method  to  determine  R.  The  counts  under  each  line  were 
summed  over  40  spectrometer  steps  (20.1*)  and  the  sums  were  corrected 
for  background;  the  background  obtained  from  the  spectral  fits  was 
used.  Then  each  sum  was  corrected  to  a  full  line  sum,  where  we 
assumed  the  profiles  were  Lorentzlan  in  the  wings.  This  is  a  good 
assumption  since  the  line  shape  in  the  wings  is  dominated  by  the  in¬ 
strumental  profile,  which  is  approxima tely  Lorentzlan.  The  correc¬ 
tions  to  the  intercomblnatlon  and  forbidden  line  sums  resulting  from 
this  process  were,  respectively,  1.011  and  1.028  times  the  resonance 
line  correction.  The  line  profiles  are  of  such  width  that  the  reso¬ 
nance  and  forbidden  lines  make  small  contributions  to  the  intercombl- 
nation  line  fluxes.  Therefore,  1.09%  of  the  corrected  resonance  line 
count  and  0*542  of  the  corrected  forbidden  line  count  were  subtracted 
from  the  corrected  intercomblnatlon  line  count.  Finally,  each  sum  was 
corrected  for  spectrometer  response  as  s  function  of  wavelength  and  R 
was  calculated.  The  results  are  shown  in  Table  1  (Rgy^).  We  also 
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show  R  from  the  spectral  fitting  program,  where  acceptable  fits  were 
obtained.  "F"  beside  the  spectrum  number  Indicates  that  the  fast 
stepping  rate  was  used. 

It  Is  difficult  to  establish  with  certainty  that  no  flare  was  in 
progress  at  a  given  time.  We  have  used  published  flare  listings  and 
measurements  by  the  MONEX  low-energy  X-ray  monitor,  also  on  P78-1,  to 
check  for  possible  flares.  There  were  reported  subflares  In  the  ob¬ 
served  regions  near  the  observation  times  for  spectra  6,  22,  and  28, 
and  possible  minor  bursts  detected  by  the  MONEX  experiment,  which 
views  the  whole  solar  disk,  for  numbers  24  and  29.  These  possible 
flares  will  not  affect  the  conclusions  below,  but,  at  least  in  the 
case  of  spectrum  28,  it  is  doubtful  that  nonflare  densities  were  meas¬ 
ured. 

-2  — 

If  the  data  are  weighted  by  cr  ,  R,  the  mean  value  of  R  from 
Table  I,  is  3.59  but  the  distribution  of  measurements  is  not  consis¬ 
tent  with  the  assumption  that  RQ  is  3.59  and  all  observed  regions  are 
in  the  low-density  limit.  Low  R  values  can  be  explained  as  arising 
from  high-density  regions,  but  high  values  of  R  cannot  be  similarly 
eliminated.  By  eliminating  all  observations  having  R  more  than 
2o  below  3.59  (spectra  10,  28),  where  a  is  the  tabulated  standard 

deviation,  we  arrived  at  E  ■  3.66,  but  the  distribution  was  still 
inconsistent  with  RQ  -  7.  Finally,  by  eliminating  all  observations 
for  which  R  +  2o  <  3.75  (spectra  5,  10,  23,  26,  28,  and  29),  we  find 
that  7  ■  3.74  and  that  the  distribution  is  consistent  with  X  ■  R0  with 
all  remaining  observations  being  in  the  low-density  limit.  Of  course. 
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higher  values  of  RQ  are  also  possible  If  some  of  Che  observations  are 

Q 

of  regions  having  densities  above  ~  3  x  103  cm  .  Actually  R0 
increases  slowly  with  increasing  temperature,  as  Blumenthal ,  Drake, 
and  Tucker  (1972)  pointed  out.  Pradhan  and  Shull  (1981)  predict  a 
slower  Increase  than  did  these  authors.  The  value  3.95  pertains  to  T 
■  Tm  -  1.87  x  10**  K  (the  temperature  of  maximum  0  VII  population), 
which  is  probably  close  to  the  average  temperature  in  the  active 
regions  we  observed.  The  densities  in  Table  1  are  derived  under  the 
assumption  that  RQ  -  3.9  ±  0.2.  Our  measurements  are  not  corrected 
for  satellite  lines.  Bhalla,  Gabriel,  and  Presnyakov  (1975)  give  the 
satellite  line  strengths  relative  to  the  resonance  line  strength, 
which  itself  Includes  unresolved  satellite  lines.  In  our  data  analy¬ 
sis  procedure,  satellite  lines  increase  the  apparent  forbidden  line 
flux  by  2.0%  of  the  resonance  line  flux  and  the  Intercombination  line 
flux  by  0.3%  of  the  resonance  line  flux  at  T  -  Tm.  For  T  >  Tm  the 
effect  is  smaller.  The  net  result  is  that  the  measured  R  is  Increased 
by  ~  1%  at  T  -  Tm. 

HI.  G,  THE  TRIPLET-TO-SINGLET  LINE  RATIO 

Acton  and  Brown  (1978)  showed  that  the  ratio  G,  defined  in  equa¬ 
tion  (1),  can  be  useful  as  an  indicator  of  the  state  of  ionization 
equilibrium.  All  three  upper  levels  Involved  in  G  are  populated  by 
recombination,  but  the  triplets  are  affected  more  than  the  singlet. 
In  a  recombining  plasma,  where  the  0  VIII  population  is  higher  than 
would  be  expected  from  the  electron  temperature,  G  is  higher  than  its 


ionization  equilibrium  value.  Conversely,  G  is  lower  in  an  ionizing 
plasma.  Unfortunately,  because  of  degradation  of  the  resonance  line 
flux  by  resonance  scattering,  a  single  measurement  of  G  is  inadequate 
to  determine  the  state  of  ionization  equilibrium  for  the  case  of  a 
recombining  plasma.  Resonance  scattering  significantly  affects  G  for 
both  0  VII  and  Ne  IX. 

Acton  (1978)  has  applied  resonance  scattering  theory  to  the  coro¬ 
nal  helium-like  ions.  The  process  consists  of  the  absorption  of  a 
photon  by  an  ion  in  the  Is*  Sq  ground  state  to  produce  an  excited  ion 
in  the  ls2p  state,  which  then  decays  to  the  ground  state.  In  our 
observations  of  active  regions,  the  following  processes  generate  or 
alte**  the  detected  X-ray  flux:  emission  in  the  active  region,  scat¬ 
tering  of  X-rays  emitted  by  the  active  region  into  the  detector,  and 
scattering  of  radiation  out  of  the  beam  headed  for  the  detector. 
Detection,  either  directly  or  after  scattering,  of  X-rays  emitted 
outside  active  regions  is  negligible.  This  is  verified  by  the  fact 
that  0  VII  and  Ne  IX  lines  are  measurable  by  the  SOLEX  B  spectrometer 
only  when  an  active  region  is  in  the  field  of  view. 

Observations  from  Sky lab  made  Jit  obvious  that  the  corona  is 
highly  structured.  Interpretations  of  the  observations  have  led  to 
the  idea  that  nearly  all  the  coronal  plasma  is  confined  in  closed 
magnetic  flux  tubes  (Rosner,  Tucker  and  Vaiana  1978;  Serio  jet_  al. 
1981).  Recently  Pallavicinl  et  al.  (1981)  found  that  the  observed 
properties  of  a  wide  variety  of  loops  within  and  Interconnecting 
active  regions  can  be  reproduced  satisfactorily  with  simple  static 
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loop  models.  The  old  model  of  a  uniform  atmosphere  in  hydrostatic 
equilibrium  is  replaced  by  an  ensemble  of  loops,  each  in  hydrostatic 
equilibrium.  In  looking  at  active  region  emission  we  are  looking 
through  a  collection  of  loops.  The  optical  depth,  t,  is  proportional 
to  path  length,  given  by  Ars(9),  where 


s(0) 


(4) 


In  equation  (4),  9  is  the  angle  between  the  line  of  sight  and  the 
inward  normal  to  the  solar  surface,  rs  is  the  solar  radius,  and  Ar  is 
the  thickness  of  the  shell  in  which  the  scattering  occurs.  When 
2Arsec^9/rs  «  1,  s(9)  *  sec9;  that  is,  a  plane  atmosphere  is  a  good 
approximation.  Equation  (4)  is  not  valid  for  0  >  90®,  and  for  0  -  90° 
the  collimator  field  of  view  is  large  enough  that  9  is  poorly  defined 
so  the  equation  should  not  be  used.  Because  of  the  nonuniform 
structure  mentioned  above,  the  relationship  between  t  and  s(0)  is 
expected  to  show  considerable  scatter. 

An  active  region  emits  resonance  line  radiation  toward  the  detec¬ 
tor  both  directly  and  by  scattering.  The  scattering  component  is 
small  because  the  optical  thickness  of  the  active  region  is  small.  We 
estimate  that  scattering  accounts  for  <_  5X  of  the  resonance  line  flux 
emanating  from  an  active  region  and  assume  that  this  fraction  does  not 
vary  with  angle.  The  estimate  is  based  on  analysis  to  be  discussed 
later  in  this  section. 
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In  order  for  an  X-ray  to  be  scattered  Into  the  detector,  the 
scattering  must  take  place  within  the  field  of  view  and  the  final 
propagation  direction  must  be  within  the  very  small  detector  solid 
angle*  Scattering  into  the  detector  becomes  less  important  as  the 

I 

angle  8  Increases.  To  see  this,  consider  a  scattering  element  of 
volume  dV  in  the  "quiet  corona"  at  a  height  z  above  an  emitting  point 
P,  such  that  the  vector  from  P  to  dV  makes  an  angle  8'  with  the  local 
vertical.  The  solid  angle  subtended  by  dV  at  P,  and  therefore  the 
flux  through  dV,  is  proportional  to  (zsecB')"^.  The  probability  of  a 
scattering  occurring  in  the  height  range  z  to  z  +  dz  is  KdzsecS', 

i 

where  k  is  the  resonance  scattering  "absorption  coefficient."  Thus 
the  number  of  scatterings  toward  the  detector  per  unit  z-range  is 
^  proportional  to  cos8*.  Since  active  regions  and  the  detector  field  of 

view  both  have  characteristic  sizes  in  the  arc  minute  range  (1  arc  min 
■  4  x  109  cm),  8'  takes  on  a  wide  range  of  values,  but  the  average 
I  value  increases  with  8.  Thus  scattering  into  the  detector  decreases 

with  increasing  8. 

If  scattering  into  the  detector  is  negligible  compared  to  scat- 
<  tering  out,  the  resonance  line  flux,  Fr,  is  proportional 

to  e  T,  where  t  is  the  optical  depth.  Then,  for  8  not  near  90°, 

•  2 - F”1  ~  eT  -  exp[s(8)  /  ic(z)dz],  (5) 

o  z 

3 

where  zg  is  the  height  of  the  source  in  the  corona,  GQ  is  the  value  G 
'  would  take  were  resonance  scattering  negligible,  and  k(z)  is  the 
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average  "absorption  coefficient"  integrated  over  the  line  profile, 
talcing  into  account  the  temperatures  of  both  the  emitting  and 
absorbing  regions.  Equation  (5)  does  not  take  into  account  changes  in 
the  line  shape  caused  by  absorption,  but  such  changes  can  be  neglected 
if  t<1,  as  it  is  here*  Taking  the  natural  logarithm  of  equation  (5) 
we  have 

m 

InG  -  [/  ie(z)dz]s(8)+  InG  .  (6) 

zg  ° 

For  8*0,  scattering  into  the  detector  may  be  comparable  in  impor¬ 
tance  to  scattering  out,  with  the  relationship  depending  on  the  geo¬ 
metry  of  the  emitting  region  and  the  detector  field  of  view.  Scatter¬ 
ing  in  becomes  less  important  with  increasing  8,  and  equation  (6)  is  a 
useful  approximation  for  determining  the  column  density  of  material 
along  the  line  of  sight. 

Table  2  gives  the  measured  values  of  G  for  both  0  VII  and  Ne 
IX.  The  spectrum  numbers  are  the  same  as  in  Table  1.  We  derived  the 
G  values  for  0  VII  by  using  the  summing  technique  in  §11.  For  Ne  IX 
the  technique  differed  in  that  each  initial  sum  was  over  30  (not  40) 
spectrometer  steps,  and  the  background  was  interpolated  from  counting 
rates  in  two  spectral  ranges  (shaded  in  Figure  2).  Figure  2  shows  the 
spectrum  around  the  Ne  IX  lines.  The  Ne  IX  intercombination  line  is 
near  the  resonance  line  and  is  much  weaker.  This  prevents  a  good  R 
determination  for  Ne  IX,  but  since  G  is  near  1.0,  corrections  to  the 
intercombination  line  are  relatively  unimportant  in  determining  G. 


TABLE  2 


G  for  O  VII  and  Ne  IX 


Spectrum 

Number 

e 

(degrees) 

s(0) 

G(0  VII) 

G(Ne  IX) 

1 

44 

1.  36 

1.  17  +  0.03 

0.78+0.02 

2 

58 

1.  75 

1.  19  +  0.03 

0.80  +  0.  03 

3 

54 

1.62 

1.  36  +  0.  04 

0.89  +  0.04 

4 

66 

2.  14 

1.82  +  0.05 

0.  89  +  0.  03 

5 

74 

2.  86 

1.  84  +  0.  04 

0.  92  +  0.  02 

6 

19 

1.  06 

1.01  +  0.02 

0.77  +  0.  05 

7 

33 

1.  18 

1.01  +0.03 

0.91  +  0.06 

8F 

19 

1.  06 

1.  08  +  0.  04 

0.81  +  0.04 

9F 

19 

1.06 

a 

0.81  +  0.03 

I  OF 

23 

1.08 

1.14  +  0.07 

0.  82  _+  0.  03 

1  IF 

32 

1.  17 

1.09  +  0.03 

0.  84  +  0.  04 

12F 

42 

1.31 

1.21  +  0.03 

0.  87  +  0.  04 

I3F 

57 

1.73 

1.00  +  0.04 

0.  86  +  0.  06 

I4F 

64 

2.  04 

1.24  +  0.04 

0.72  +  0.  05 

15F 

70 

2.47 

1.54  +  0.09 

0.90  +  0.  09 

I6F 

82 

3.  94 

1.78  +  0.08 

l.  09  +  0.  11 

17 

92 

1.93  +  0.03 

0.  98  +  0.  05 

18 

58 

1.75 

1.47  +  0.  03 

0.91  +0.05 

19 

68 

2.  30 

1.60  +  0.  03 

a 

20 

92 

2.25  +  0.04 

0.99  +  0.04 

21 

101 

2.01  +  0.04 

0.  90  +  0.  05 

22 

57 

1.72 

1.45  +0.03 

1.  07  +  0.  04 

23 

74 

2.  87 

2.25  +  0.05 

1.  01  +  0.  04 

TABLE  2  (Con't) 

G  for  O  VII  and  Ne  IX 
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Spectrum 

Number 

e 

(degrees) 

s  (0) 

G(0  VII) 

G(Ne  IX) 

24 

48 

1.  46 

1.21  +  0.  03 

0.82  +  0.04 

25 

33 

1.  18 

1.27  +  0.  03 

0.77  +  0.  04 

1 

26 

31 

1.  15 

1.  20  _+  0.  03 

0.65  +  0.06 

27 

38 

1.25 

1.  16  +  0.  04 

0.  85  +  0.  08 

28 

64 

2.  05 

1.43  +  0.03 

b 

( 

29 

71 

2.54 

1.65  +  0.03 

1.  00  +  0.  03 

30 

77 

3.20 

1. 87  +  0.  04 

0.  92  +  0.  03 

^Data  not  usable. 

°  Probable  flare,  possible  Fe  XIX  emission. 
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X-ray  Spectrum  of  an  Active  Region  Showing  Ne  7.X  Emission 
Lines.  The  regions  in  which  background  was  measured  are 
shaded. 


For  Ne  IX  such  corrections  lower  G  by  0.059,  and  for  0  VII  G  is 
lowered  by  0.015.  Errors  in  G  associated  with  the  line  shape  will  be 
smaller  than  these  corrections.  The  G  values  are  not  corrected  for 
satellite  lines.  It  is  customary  to  include  blended  satellite  lines 
in  theoretical  calculations  of  G  (Acton  and  Brown  1978;  Pradhan  and 
Shull  1981).  The  calculations  therefore  correspond  to  the  line  ratio 
that  is  actually  observed,  and  direct  comparisons  between  theory  and 
observation  are  facilitated.  For  Ne  IX  there  is  also  the  possibility 
of  a  blend  with  Fe  XIX  lines.  The  Fe  XIX  population  peaks 
at  T  ■  7  x  106  K  and  is  small  below  -  4  x  1q6  k  (Jacobs  jst_  al. 
1977).  Therefore,  we  expect  Fe  XIX  emission  from  nonflaring  plasmas 
to  be  weak.  We  have  not  found  any  evidence  of  its  existence.  The 
strongest  Fe  XIX  lines  are  blended  with  the  weak  Ne  IX  lntercomblna- 
tion  line  in  our  flare  spectra  (Bromage  and  Fawcett  1977;  McKenzie  et 
al.  1980b),  so  Fe  XIX  emission  should  be  manifested  by  a  decrease  in 
the  ratio  of  forbidden  to  intercombination  line  counts.  Thus  if  the 
measured  G  has  been  Increased  by  Fe  XIX  emission,  G  should  be  nega¬ 
tively  correlated  with  the  above  ratio.  We  find  no  such  correla¬ 
tion.  Table  2  shows  significant  variations  in  G  for  both  0  VII  and  Ne 
IX.  The  quoted  errors  are  due  to  counting  statistics  only. 

In  Figures  3  and  4  we  plot  InG  vs  s(9)  for  0  VII  and  Ne  IX, 
respectively,  omitting  points  for  which  9  £  90°.  Table  2  shows  that 
the  omitted  points  have  large  t,  as  expected.  Equation  (4)  shows  that 

g 

s(9)  depends  on  an  assumed  value  for  Ar.  We  assume  Ar  ■  4.5  x  10  cm, 
which  is  equal  to  the  scale  height,  kT/mHg9,  where  gQ  is  the. 
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Figure  3.  Plot  of  InG  vs.  s(0)  for  0  VII  Observations  of  Solar 

Active  Regions.  The  optical  depth  for  resonance  scatter¬ 
ing  is  proportional  to  s(8)  (see  text).  The  scale  at  the 
bottom  allows  G  to  be  related  to  6,  the  angle  between  the 
line  of  sight  and  the  Inward  normal  to  the  solar  surface. 
For  8  <  50.2°,  s(fl )  •  sec  6. 
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acceleration  of  gravity  at  the  Sun's  surface,  Dy  is  the  proton  mass, 
and  T  -  1.5  x  10^  K.  The  results  of  the  analysis  below  depend  weakly 
on  Ar.  The  best  fits  to  equation  (6)  are  also  shown  in  the  figures. 
Considerable  scatter,  arising  from  the  inhomogeneity  of  the  corona,  is 
evident. 

The  curves  in  Figures  3  and  4  have  slope 

m 

b  •  /  «(z)dz  .  (7) 

z 

s 

Acton  (1978)  gives  the  following  expression  for  ico(z)  the  line  center 
absorption  coefficient: 


■c  (z) 
o 


9.31  x  10 


-18 


f,  A. 
lu  i 


ne(z)l 


/  M  \1/2  -1 

(r^yj  “  • 


(8) 


In  equation  (7)  flu  is  the  absorption  oscillator  strength,  Ag  is  the 
element  abundance  relative  to  hydrogen,  A^  is  the  fractional  popula¬ 
tion  of  0  VII  or  Ne  IX,  \  is  the  wavelength  in  A,  M  the  ion  mass  in 
atomic  mass  units  and  Tg  the  temperature  in  units  of  10^  X.  The 
effective  absorption  coefficient,  derived  by  taking  the  line  shapes 
into  account,  is  nKQ(z),  where 


a 


(9) 


4 


c 

3 


4 


4 


4 


and  Te  and  T^  are  Che  temperatures  in  the  emitting  and  absorbing 
regions,  respectively. 

The  curves  can  be  used  to  evaluate  the  electron  column  densities 
in  the  quiet  corona  and  a  "typical"  active  region.  We  assume  that  the 
quiet  coronal  temperature,  Tq,  is  1.5  x  106  K  (McKenzie  £t  al.  1978; 
if  Tq  *  1.0  x  10®  K  the  results  below  would  be  changed  very  little) 
and  that  the  active  region  temperature  is  2.5  x  10®  K.  A^e  is  assumed 
to  be  6  x  10-5  and  Aq  is  taken  to  be  3  x  10~^.  These  values  are  based 
on  Acton's  (1978)  recommendation  that  Withbroe's  (1976)  0  abundance  be 
multiplied  by  0.65  and  on  the  Aq/AN6  ratio  of  Acton,  Catura,  and  Joki 
(1975).  Taking  bg  ■  0.290  and  bjje  ■  0.090  from  Figures  3  and  A,  re¬ 
spectively,  we  can  solve  for  the  electron  column  densities,  Nq  and  NA, 
in  the  quiet  corona  and  the  typical  active  region  (for  0  -  0  ).  We 
find  Nq  •  3  x  10*®  cm-^  and  NA  •  5  x  10*®  cm-^.  The  total  column 
density  through  the  active  region  is  ~  2NA«  With  a  typical  loop 
height  of  2  x  10^  cm  (Pallaviclni  et_  al.  1981)  the  average  active 
region  density  is  -  5  x  107  cm  .  If  we  estimate  the  quiet  coronal 
path  length  as  the  scale  height,  kT/myg^,  we  have  ne  at  the  base  of 
the  quiet  corona  ~  7  x  10®  cm“®.  Thus  the  active  region  density  is 
approximately  an  order  of  magnitude  higher  than  that  in  the  quiet 
corona.  With  these  values,  less  than  ~  5%  of  the  0  VII  emission  and 
less  than  ~  0.5%  of  the  Ne  IX  emission  arise  from  outside  the  active 
region.  Approximately  80%  of  the  0  VII  absorption  and  ~  40%  of  the  Ne 
IX  absorption  take  place  in  the  quiet  corona. 
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If  scattering  into  the  detector  could  be  neglected,  GQ  could  be 
found  by  extrapolating  the  curves  in  Figures  3  and  4  to  s(8)“  0  (i. 
e. ,  to  T  ■  0).  However,  scattering  into  the  detector  is  not  negli¬ 
gible.  Its  importance  relative  to  scattering  out  depends  on  geo¬ 
metrical  relationships  between  the  emitting  material  and  the  instru¬ 
ment  field  of  view.  Scattering  is  least  important  for  6*0.  Judging 
from  sample  calculations,  raster  maps  that  give  us  geometrical 
Information,  and  G  values  for  small  0,  we  estimate  that  GQ  (0  VII)  • 
1.0  ±  0.1.  For  Ne  IX,  G0  can  be  estimated  from  G(6*  0)  since  Ne  IX 
scattering  is  almost  negligible  when  8  is  small.  An  average  of  data 
from  spectra  6-11  and  25-26  with  data  weighted  by  o  (from  Table  2) 
gives  G0(Ne  IX)  ■  0.80  ±  0.05.  The  estimated  GQ  (0  VII)  agrees  with 
Pradhan  and  Shull's  (1981)  calculated  value  of  1.00,  but  our  GQ(Ne  IX) 
differs  substantially  from  their  value  of  1.05. 


IV.  DISCUSSION 

The  0  VII  line  ratio,  R,  can  be  used,  as  in  equation  (2),  to 
measure  electron  densities  in  the  solar  corona  at  T  *  2  x  10^  K. 
Uncertainties  in  ne  arise  from  uncertainties  in  R  and  in  RQ.  Errors 
in  Rq  are  particularly  Important  near  the  threshold  density  for  which 
the  ratio  is  applicable.  Dividing  equation  (2)  by  equation  (3)  gives, 
for  the  uncertainty  arising  from  ARq  alone, 


An  AR 
e  o 

—  “i r=T  » 

e  o 
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where  A  denotes  “uncertainty  in."  In  flares,  usually,  ARq<  <  RQ  -  R 
(McKenzie  et  al.  1980a;  Doschek  et  al.  1981),  so  the  uncertainty  in  RQ 
is  unimportant.  In  nonflaring  spectra,  ARq  is  frequently  comparable 
to  RQ  -  R,  and  uncertainties  in  ne  are  consequently  large.  A  major 
aim  of  the  study  described  here  has  been  to  evaluate  RQ. 

From  the  analysis  of  29  0  VII  spectra,  we  find  that  the  data  are 
consistent  with  RQ  >  3.74.  Correction  for  blended  satellite  lines 
would  reduce  this  lower  limit  by  ~  1%.  We  have  implicitly  assumed 
that,  for  our  measurements,  T  -  Ta*  It  is  likely  that  there  was  a 
small  variation  in  temperature  among  the  various  observed  regions.  We 
estimate  that  RQ  varied  by  ~  32  (la).  The  R0  range  we  arrive  at  can 
be  considered  to  be  applicable  to  nonflare  spectra.  During  flares  RQ 
might  be  higher,  but  densities  are  usually  high  enough  that  uncertain¬ 
ties  in  RQ  do  not  result  in  major  errors  in  n0.  We  recommend  that  the 

same  value  of  R0  be  used  in  analyzing  all  0  VII  data. 

In  making  the  density  determinations  in  Table  1  we  assumed  Rq  * 

3.9  i  0.2.  Uncertainties  in  R  and  in  R0  were  taken  into  account.  Our 
data  allow  R0  as  low  as  3.7,  and  all  of  the  data  that  we  know  of  are 
consistent  with  RQ  _<  4.1.  Although  the  uncertainty  in  RQ  is 
small,  ~  52,  the  uncertainties  in  ne  in  Table  1  are  large.  In  only 
one  case,  spectrum  28,  is  An  /n  <  0.4,  and  in  this  case  the  region 
observed  may  have  been  flaring.  However,  even  a  density  determination 
with  ~  502  uncertainty  may  be  valuable  in  studying  the  physics  of 
coronal  active  regions. 
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The  attenuation  of  resonance  line  flux  by  resonance  scattering  is 
apparent  for  both  0  VII  and  Ne  IX.  In  fact  it  appears  that  the  0  VII 
flux  is  almost  always  subject  to  significant  attenuation.  The  reso¬ 
nance  line  flux  alone  has  limited  usefulness  even  in  deducing  emission 
measure.  Uhere  resonance  line  flux  alone  is  available,  it  should  be 
corrected  by  multiplying  it  by  G/GQ,  where  G  is  determined  from  Figure 
3  as  a  function  of  8  and  GQ  -  1,0.  For  Ne  IX,  resonance  scattering  in 
the  emitting  region  is  as  Important  as  in  the  Intervening  corona. 
This  makes  it  more  difficult  to  correct  the  Ne  IX  flux  for  resonance 
absorption.  On  the  other  hand,  since  resonance  scattering  is 
relatively  weak  for  Ne  IX,  the  use  of  the  curve  in  Figure  4  in 
correcting  the  resonance  line  flux  will  usually  limit  errors  to 
<  201.  Such  errors  are  not  serious,  given  the  current  state  of 
absolute  X-ray  spectrometry  of  the  solar  corona. 

Our  Interest  in  G  was  motivated  by  its  potential  usefulness  as  an 
Indicator  of  the  state  of  ionization  equilibrium  of  the  emitting 
plasma.  Because  of  resonance  scattering,  a  single  measurement  of  G  is 
Inadequate  to  determine  whether  or  not  the  plasma  is  in  ionization 
equilibrium.  However,  for  0  VII  the  resonance  scattering  takes  place 
primarily  in  the  corona  between  the  source  and  the  observer.  This  is 
probably  especially  true  for  flares,  where  the  temperature  is  high  and 
the  0  VII  population  is  relatively  low.  If  a  preflare  determination 
of  G  is  available  and  the  flare  values  are  unchanged,  this  can  be 
taken  as  evidence  that  large  departures  from  ionization  equilibrium 
are  not  occurring.  We  have  analyzed  only  a  few  flares  for  which 
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preflare  G's  are  available,  but  in  no  case  have  we  observed  a 
significant  variation  in  G  during  a  flare.  This  result  is  in  accord 
with  the  conclusion  of  Doschek,  Feldman,  and  Cowan  (1981)  that  the 
high-temperature  plasmas,  from  which  Fe  XXII  -  Fe  XXV  lines  are 
emitted,  are  in  ionization  equilibrium. 
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laboratory  operations 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting  exper¬ 
imental  and  theoretical  investigations  necessary  for  the  evaluation  and  applica¬ 
tion  of  scientific  advances  to  new  military  space  systems*  Versatility  and 
flexibility  have  been  developed  to  a  high  degree  by  the  laboratory  personnel  in 
dealing  with  the  many  problems  encountered  in  the  nation's  rapidly  developing 
space  systems.  Expertise  in  the  latest  scientific  developments  is  vital  to  the 
accomplishment  of  tasks  related  to  these  problems.  The  laboratories  that  con¬ 
tribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  vehicle  and  reentry  aerodynamics  and  heat 
transfer,  propulsion  chemistry  and  fluid  mechanics,  structural  mechanics,  flight 
dynamics;  high-temperature  thermomechanics,  gas  kinetics  nd  radiation;  research 
in  environmental  chemistry  and  contamination;  cw  and  pulsed  chemical  laser 
development  including  chemical  kinetics,  spectroscopy,  optical  resonators  and 
beam  pointing,  atmospheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmo¬ 
spheric  optics,  light  scattering,  state-6pecif ic  chemical  reactions  and  radia¬ 
tion  transport  in  rocket  plumes,  applied  laser  spectroscopy ,  laser  chemistry, 
battery  electrochemistry,  space  vacuum  and  radiation  effects  on  materials,  lu¬ 
brication  and  surface  phenomena,  thermionic  emission,  photosensitive  materials 
and  detectors,  atomic  frequency  standards,  and  bloenvlronmental  research  and 
monitoring. 

Electronics  Research  Laboratory:  Microelectronics,  CaAs  low-noise  and 
power  devices,  semiconductor  lasers,  electromagnetic  and  optical  propagation 
phenomena,  quantum  electronics,  laser  communications,  Ildar,  and  electro-optics; 
communication  sciences,  applied  electronics,  semiconductor  crystal  and  device 
phy*<_s,  radiometric  imaging;  mi llioeter-wave  and  microwave  technology* 

Information  Sciences  Research  Office:  Program  verification,  program  trans¬ 
lation,  performance-sensitive  system  design,  distributed  architectures  for 
spaceborne  computers,  fault-tolerant  computer  systems,  artificial  intelligence, 
and  microelectronics  appl ications. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metal  matrix 
composites,  polymers,  and  new  forms  of  carbon;  component  failure  analysis  and 
reliability;  fracture  mechanics  and  stress  corrosion;  evaluation  of  materials  In 
space  environment;  materials  performance  in  space  transportation  systems;  anal¬ 
ysis  of  systems  vulnerability  and  survivability  In  enemy-induced  environments. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radiation 
from  the  atmosphere,  density  and  composition  of  the  upper  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation  of 
plasma  waves  in  the  magnetosphere;  solar  physics,  infrared  astronomy;  the 
effects  of  nuclear  explosions,  magnetic  storms,  and  solar  activity  on  the 
earth’s  atmosphere,  ionosphere,  and  magnetosphere;  the  effects  of  optical, 
electromagnetic,  and  particulate  radiations  in  space  on  space  systems. 


